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Abstract 

We develop a method of computing the excited state energies in Integrable Quantum 
Field Theories (IQFT) in finite geometry, with spatial coordinate compactified on a cir- 
cle of circumference R. The IQFT "commuting transfer-matrices" introduced in 0] for 
Conformal Field Theories (CFT) are generalized to non-conformal IQFT obtained by per- 
turbing CFT with the operator $1,3. We study the models in which the fusion relations for 
these "transfer-matrices" truncate and provide closed integral equations which generalize 
the equations of Thermodynamic Bethe Ansatz to excited states. The explicit calculations 
are done for the first excited state in the "Scaling Lee- Yang Model" . 
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1. Introduction 

Much progress in understanding 2D Integrable Quantum Field Theories (IQFT) has 
been achieved during the last two decades (see e.g. for a review). In particular, many 
models have been exactly solved "on-shell" i.e. their mass spectra and their factorizable 
5'-matrices have been determined. Important progress was also made towards "off-shell" 
solution of IQFT. The approach based on so called Thermodynamic Bethe Ansatz (TBA) 
technique has proven to be very successful . This approach allows one to calculate the 
ground state energy Eq[R) of IQFT in the finite-size geometry, where the spatial coordinate 
is compactified on the circle of circumference R, provided the "on-shell" solution is known. 
The problem being reduced to solving certain non-linear integral equation (TBA equation). 
The quantity Eq{R) makes it possible to probe the scale dependence in the IQFT, in 
particular, to find in the limit — > the central charge c of the Conformal Field Theory 
(CFT) which governs the short distance behavior of the IQFT. 

So far the power of TBA approach was limited to the ground state energy Eq{R) 
(or the ground state energies in "twisted" sectors 0, [0). At the same time it seems 
to be interesting to calculate the full energy spectrum En{R) of the finite-size IQFT. 
These energies would provide interpolation between the CFT spectrum (determined by 
the CFT central charge c and the conformal dimensions A^) at i? and massive 
spectrum (determined by the particle masses and the S'-matrix amplitudes) at i? ^ oo. 
However it is not clear how the ideas behind traditional derivation of TBA equation could 
be generalized to incorporate the excited states. 

In our recent paper [|l|] we have shown how one can define the family of operators 
Tj(A); j = 1/2,1,3/2,... — the "QFT transfer-matrices" S which act in the space of 
states of the finite size IQFT and commute for different values of the parameter A 

[T,(A),T,.(A')] = 0. (1.1) 

We have also shown that the "fusion relations" for the operators Tj(A) 

T,{qH)T,{q--^X) = I + T,_i/2(A)T,+i/2(A) . (1.2) 

^ These operators appear as the continuous QFT versions of Baxter's commuting transfer- 
matrices 0, and therefore we maintain using this term although the original meaning of the 
term "transfer-matrix" here is apparently lost. 
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(here g is a parameter related to the Virasoro central charge c), give rise to the functional 
equations for the eigenvalues of these "transfer-matrices" and that for the case of the 
ground state eigenvalue these reduce to the TBA equations; at the same time the func- 
tional equations apply to the excited states as well and following the method developed by 
Kliimper and Pearce 0] one can derive from them the nonlinear integral equations which 
determine in principle the excited states En{R)- 

Our analysis in [|l| explicitly concerns the case of CFT (more precisely, the c < 1 
"minimal" CFT [0), where the space of states is decomposed in the sum of the direct 
products of "right" and "left" irreducible highest-weight Virasoro module Va„ and the 
operators T(A) act in each of these spaces Va^ separately. However, one can consider 
more general non-conformal IQFT which are obtained by perturbing the "minimal CFT" 
with the operator $1^3 [|ll|. In this paper we show that the commuting operators similar 



to Tj(A) can be constructed in the perturbed theory as well (we denote these "perturbed" 
transfer matrices as Tj(/i|A), where is the parameter of the perturbation), and that the 
most important properties of the operators Tj{X) obtained in |[^], including the fusion 
relation ( p. .21) , remain valid for the "perturbed" operators Tj(|u|A). Although the analytic 
properties of the operators Tj{^\X) as the functions of A turn out to be significantly different 
from those of the "conformal" operators Tj(A), it is still possible to derive the nonlinear 
integral equations which determine the finite-size energy levels E^iR) in the perturbed 
theory. 

To demonstrate the efficiency of our approach we calculate the energy Ei{R) of the 
first excited state of the so called Scaling Lee- Yang Model (SLYM). This model is obtained 
by perturbing the c = —22/5 "minimal" CFT A^2/5 with the operator $ — $1^3 of the 
dimension A = —1/5, i.e. it is defined by the action 

AsLYM = •^M2/5 + '"^ y ^^^^ ' ^^-^^ 

where AM2/5 action of the CFT M.2/5 A is the coupling parameter. The QFT 



defined by (|1.3|) is integrable |]Tl| and its on-shell solution (i.e. the particle content and the 



factorizable S-matrix) was found in |]T2|. The theory contains one sort of scalar particles 
with the mass m and their two-particle S-matrix is 

^ sinh(g) + i sin(7r/3) 
^ ' sinh(0) -zsin(7r/3) ' ^ ' 
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where 9 = 9i —6*2; Oi and 6*2 are the rapidities of the scattering particles (i.e. their momenta 
are pi = m sinh(6'i), p2 = m sinh(6'2) )• The finite size ground state energy Eq{R) of ( pT^ ) 
was obtained in using TBA technique. Also, few first excited levels En{R) was studied 
on the basis of "Truncated Conformal Space method" in |T^. More recently the exact 
relation between the mass m and the coupling parameter jl in ( |1.3| ) was found in [|n 

.,_. 2V^53/4 (r(2/3)r(5/6))^^/^ 

^ 167r6/5 r(3/5)r(4/5) ' ^ ^ ^ 

We apply the approach based on the functional equations for the transfer-matrices T to 
calculate the finite-size energy of the first excited state, the state which is interpreted as 
the one-particle state (with the zero-momentum particle) at — 00 and approaches the 
CFT state associated with the identity operator of the CFT as R 0. 

The paper is organized as follows. In Sect. 2 we summarize the properties of the 
operators Tj(A) of the Conformal Theory defined in |jl|]. The "transfer-matrices" Tj{n\X) 
of the perturbed theory are introduced in Sect. 3 where also our basic conjectures about 
the analytic properties of these operators and their relation to the local Integrals of Motion 
(IM) are formulated. In Sect. 4 we study the simplest eigenvalues of the operator Ti/2(-^) 
in the CFT A^2/5- Similar analysis of the eigenvalues of Ti/2(a*I-^) SLYM is carried out 
in Sect. 5 where we derive the excited state energy Ei{R). In Sect. 6 the relation to the 
sine- Gordon theory is discussed and further applications of our approach are outlined. 



2. The operators Tj{X) in CFT 

In this section we briefiy summarize the properties of the operators Tj(A) in CFT as 
they are defined in |^. 

The space of states in a CFT decomposes as 

HcFT = ©a(VA. ®VaJ , (2.1) 

(for simplicity we consider here the "diagonal" CFT) where Va^ is the irreducible highest 
weight representation of the "left" Virasoro algebra Vir 

C n 

[L„, Lm] = (n-m) L^+m + ^2 ~ ^n+m,0 , (2.2) 

with the highest weight A^. Correspondingly, Va^ stands for the representation space of 
the "right" Virasoro algebra Vir identical to Vir. The number c in ( p.2|) is known as the 
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"central charge". The "left" Virasoro generators are expressed in terms of the "left" 
chiral component of the energy-momentum tensor T{z) {z = x + iy) as 



c . R dx 



Here we assume that the CFT is defined in the geometry of an Euclidean cylinder so that 
y is interpreted as the Euclidean time and the spatial coordinate x is compactified on a 
circle of the circumference R. The field T{z) therefore satisfies the periodicity condition 
T[z + i?) = T[z). The "right" Virasoro generators Ln are defined in terms of the "right" 
chiral component T{z), z = x — iy hj similar integrals. 



Any CFT possesses infinitely many local Integrals of Motion (IM) l2fc-i ||15|| , [p!6| 
which can be written as 



I 



2fc-l 



dx 

/ —T2kix + iy), (2.4) 
Jo 27r 



where T2k{z) are certain local fields, polynomials in T{z) and its derivatives. For example 

T2{z) = T{z), T,{z) =: T\z) :, T,{z) =: T\z) : +^ : {T'{z)f 

T2k{z) =: T^{z) : + terms with the derivatives , 

Here : : denote appropriately regularized operator products, see [Q for details. There 
exists infinitely many densities (^]^) (one for each integer k [^) such that all the integrals 
( P^) commute 

[hk-uhi-i] = . (2.6) 

The "transfer-matrices" Tj are most conveniently defined in terms of Feigin-Fuchs 
free field representation of the Virasoro algebra []T^, Introduce the "left" chiral Bose 
field! 

<f>iz)=Q+^-^P-^Yl ^e^"^-^ , (2.7) 
it ^ n 

n=lO 

where the operators Q,P,an satisfy the canonical commutation relations 

? 71 

[Q, ^] = 2 [^-' =2^^ '^"+"^'° ^^-^^ 

and /3 is a parameter which will be defined below. Also, let Tp be the Fock space, i.e the 
space generated by the action of with n < on the "vacuum" state | p) which satisfies 
the equations 

P \ p) = p \ p) , On I p) = for n > . (2.9) 



This field differs in normalization from tlie field if which is used in Q, ip{z) 
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Note that ( p77| ) implies that the field (f){z) is quasi-periodic, (f){z + R) = (j){z) + 2tx P. As is 
well known the energy-momentum tensor 

T{z) = : {d.4>? : +z (1 - Z?"^) dlcj) - ^ (2.10) 



generates the Virasoro algebra ( |2.2| ) with the central charge 

c= 13-6(/32-K/?-2) , (2.11) 

which therefore is represented in the Fock space Tp. We will always assume that c < 1 
and hence that 13 is real. For generic values of /3 and p the space JFp is isomorphic to the 
irreducible highest weight module Va with 

while for particular values of these parameters, when the "null-vectors" appear (so called 
"degenerate representations" ) , the irreducible representation Va can be obtained from jFp 
by factoring out the corresponding invariant subspaces. In fact, in what follows we always 
will imply that all such an invariant subspaces in jFp are factored out and Vp will stand 
for the "Virasoro irreducible Fock space" , Vp = JFp / (invariant subspaces) . With this 
convention the space Vp is isomorphic to Va(p)- 

Let j be a non-negative integer or half-integer number. Consider the following (2j -\- 
1) X (2j -fl) matrix 

Lj(A) = e^^^^^ Pexp dz{: e'^'^^^'^ : Ej+ : e^'"^^^^ : q-^F^)^ , (2.13) 

where Ej, Fj and Hj stand for the 2j + 1 dimensional representation nj of the quantum 
algebra t/g(s/(2)) with 

q = e*"^' , (2.14) 



and the symbol V in ( 2.13 ) denotes the "path ordering" of the both operator and matrix 
products. The matrix elements [Lj(A)]^, A,B = — j, —j + 1, j are operators which act 
in the space 

= ®u=-oo^P+nl3^ ■ (2.15) 

Note that the exponential fields : e^^*'^*^^-' : in ( |2.13|) have anomalous dimensions and 
therefore the parameter A carries the dimension 

A ~ [lengthY^~^ . (2.16) 
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As discussed in the expression (|2.13| ) applies directly only for < 1/2 (i.e. for 

c < -2); for [5"^ > 1/2 the integral in (|2J3D can be defined through appropriate analytic 
continuation in 

The operator matrices Lj(A) satisfy the Yang-Baxter relation 



J2 [^..'(A/A')]$$; [L,(A)]^ [L,.(A 

C=-j C'=-j' 



c 



C=-j C'=-j' 



CC ' 



where [-Rji'(A)]^^, are the elements of the i2-matrix acting in the product of 2j + 1 and 
2/ + 1 dimensional representations of t/q(sZ(2)). As a simple consequence of ( p.l7|) the 
"transfer- matrices" 



(2.17) 



T,(A)=tr.Je^-^^^ L,(A)] = J2 [e^'^^^^ L,(A) ] 

A=-j 



(2.18) 



form a commuting family of operators, i.e. they satisfy (1J-). Note that although the 
off-diagonal elements of the matrices ( |2.13|) act between different components of (|2.15| ) the 
operators ( p.l8|) invariant ly act in the spaces Vp, 



[2.19) 



It is possible to show using (|2.13|) , (|2.17|) that the operators Tj (A) satisfy the fusion 
relations (|1.2|) [|l|, which allow one to express any operator Tj(A) with j > 1/2 through 
the "fundamental" one T(A) = Ti/2(A). 

The operator T(A) is an entire function of A^ (the same is true for all the operators 
Tj(A)). Its power series expansion in the variable A^ follows directly from the definitions 

(ra,(ra 

oo 

T(A) = 2 cos(27rP) + ^ A^'^ Gk , (2.20) 



k=i 



where Gk are basic "nonlocal IM" 



Gk=q' 



X>H(xi,...,X2fc) 

_j_ g-2i7rP . g2i0(xi) .. g-2i<^(x2) .. g2i<^(x3) . . g-2j0(x2fe) . 



(2.21) 



Here the integration is performed over the domain Dfi[xi, ...,X2k) = {R > xi > X2 > •■• > 
a^2fc > 0}. In the opposite hmit ^ oo the operator T(A) admits asymptotic expansion 

i, m 

oo 

logT(A) - {KR/2n) X^^+^^ I - ^ A(i-2-)(i+?) h^^-i . (2.22) 



n=l 



in powers of A (^+^) with the coefficients proportional to the local IM ( |2.4[ ). The numerical 
factors K and are given by 



2v^ - I) / , ox 

r(i-l) V ^ ^ ^ 



(2.23) 



ni 



Here and below 



r(i + (n-i)0 



The above expansions suggest that the eigenvalues of the operator T(A) unite all the 
information about the eigenvalues of both nonlocal IM ( p.21|) and local IM (|2.4|) in a very 
interesting manner. In what follows we consider almost exclusively the vacuum eigenvalues 



T{X)\p)=T{t,p)\p) 



t = y 



2n 



(2.25) 



According to ( |2.2C1| ) and ( |2.22| ) this function expands both in convergent series around 
t = and in asymptotic series in the vicinity of t = oo, 



T{t,p) = 2cos{2'Kp) + Y.^^ Grip) 



k=i 



oo 

exp (/^ Cn f(V2-n)(l+0 ^(A(p))) 



(2.26o) 



(2.266) 



n=l 



where G^°'^{p) and /2^1]^(A) are, respectively, the vacuum eigenvalues of the nonlocal IM 
( ^^211) and the local IM (U) with R = 2ti. Note, that the expression in (|2.266|) is, in fact, 
an asymptotic series with zero radius of convergence. The eigenvalues G'"^"'^{p) are given 
by the integrals 



n 

n 



2tt 



Ul 



dui / dvi 



Vl 



U2 



dvn 



dUr 



dVr 



J>1 
n 



2 sin( 



2/3-' 



) 2sin( 



Vi - V4 



2/3-' 



n(2^'"(^) 

j>i 



-2(3^ 



X 



n(2sin(^^)) ' 2cos(2p(. + 5^( 



-20 



j>i 



1=1 



(2.27) 
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while the eigenvalues ^^-li^) can be calculated (in principle) from the explicit formulas 
for the associated local densities T2n{z) in ( |2.4D . Unfortunately for general k neither 
the integrals (|2.27| ) can be calculated explicitly nor the densities T2k{z) are known in a 



closed form. Some partial results about the vacuum eigenvalues G'^"'^(p) and l2n-i{^) are 
collected in Appendices A and B. 

For 0^ < 1/2 the entire function T(t,p) is completely determined by the positions tk 
of its zeroes in the complex t-plane. The locations of these zeroes are known exactly in 
two particular cases. First, in the "classical limit" 0, (with p and A kept fixed) the 

explicit formula for the vacuum eigenvalue T{t, p) can be easily obtained directly from the 
definition (|T|)i 

T(t, p) = 2 cos (27rv/p2 - t ) ; /^^ = . (2.28) 



In this case 



(2k + l)"^ 

tk=P^- ^ , A; = 0,1,2,.... (2.29) 

lb 



Another case where T{p,t) is known explicitly is the "free-fermion point" /3 = 1/2 [|2^ 



2tt e^^^ 

r(l/2 -\-2p + Tit) r(l/2 — 2p + Txt) 
where C is some constant. This function has two sequences of zeroes 

) = i (±2p - 1/2 - /c) , /c = 0, 1, 2, ... . (2.31) 

For generic < < i the zeroes tk do not have any simple form. Nonetheless a brief 
inspection of ( p.26a| ) and ( |2.29| ) , ( |2.31| ) makes it natural to adopt the following 

Conjecture 1. For real p and < < 1/2 all zeroes tk ofT(t,p) are real. There 
are infinitely many negative zeroes which accumulate towards t = —oo; for given value p 
there are also exactly Int(2|p| + 1/2) positive zeroes (lnt{x) denotes integer part of x). In 
particular, if \p\ < 1/4 all zeroes tk are negative. 

This assumption will be used in Sect. 4 in deriving closed nonlinear integral equation 
for the vacuum eigenvalues in particular case — 2/5. 



^ To derive (2.28) one observes that the oscillator modes can be ignored in the vertex oper- 
ators : e^^*'^ : as they do not contribute to the vacuum eigenvalues in this limit. Also, in this limit 
the zero mode operators P and Q commute and one obtains Tj{t,p) = tvj ^e^'" iVtEj+VtFj+ipHj) 
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3. The Operators Tj{^\X) in the $13 Perturbed Theory 

Our previous discussion was concerned with the Conformal Field Theory. It is well 
known that one can obtain more general non- conformal Integrable Quantum Field Theory 
by perturbing certain CFT with appropriately chosen relevant local operator. The most 
known example is a "minimal CFT" perturbed with the degenerate field $1,3, i.e. the 
primary field with the conformal dimensions 

Ai,3 = Ai,3 = -1 + 2/32 = 1 - ^ , (3.1) 



which satisfy the third-level null vector equations [10 



(2 (1 + 2/32) L_3 - 4 L_iL_2 + Ll,) $i,3(^, = , 
(2 (1 + 2/^2) L_3 - 4 L_iL_2 + H,) $i,3(^, z) = . 



(3.2) 



We will assume the canonical normalization of this field 

{c^^,s{z,z)^^,s{z',z'))^[{z-z'){z-z')]-''^''' as {z, z) ^ {z' , z') . (3.3) 
The perturbed theory is defined by the formal action 

ApcFT = AcFT + y d'^x $1,3(2:, z) , (3.4) 

where Ac ft denotes the "action" of CFT and jl is the coupling parameter with the 
dimension 

fi ~ [length]"^^'''^'^ = [length]^^' ~^ . (3.5) 

As in Sect. 2 we assume that the QFT is defined in the geometry of a cylinder with the 
spatial coordinate x = z compactified on a circle with the circumference R. Correspond- 
ingly, the field $1,3 in ( |3.4|) is taken to be periodic function of x, i.e. $1,3(2 + R, z + R) = 
$1,3(^5 z). As was shown in [11| the QFT defined by (3j^) is integrable in the sense that 



it possesses infinitely many local IM which appear as certain deformations of the local IM 



of the original CFT. Although for fi ^ the fields T2k (p.5|) are no longer holomorphic 



fields (i.e. they depend both on z and z) they satisfy the continuity condition 

dzT2k{z, z) = fi^ dz<c>2k-2iz, z) , (3.6) 
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where 02k-2{z, z) are certain local fields. The right-moving chiral sector of the original 
CFT contains also the antiholomorphic fields T2k', in the perturbed theory (|3.4| ) they satisfy 
similar continuity equations 



d,T2k{z, z) = fi' d-,Q2k-2{z, z) . (3.7) 



It follows that the integrals 



f dx 

hk-i = — {T2k{x + iy,x- iy) + /i^ Q2k-2{x + iy,x- iy)) , 
^2k-i= / — iT2k{x + iy,x -ly) fi^ Q2k-2{x + iy,x -ly)) 

Jo 

in fact do not depend on the "Euclidean time" y, i.e. they are Integrals of Motion (IM). 



It is possible to show that the operators (q^) commute among themselves 



[l2fc-l,l2Z-l] = [hk-lM-l] = [hk-iM-i] = . (3.9) 

Note that the field /i^ Oq = /i^ Oq = (1 — ^1,3) ^1,3 coincides with the trace of the 
energy-momentum tensor T^^p of the QFT (|3.4| ), and the IM Ii and Ii are related to the 
Hamiltonian EI and the spatial momentum P of this QFT 

e = ii p = ii -Ii . (3.10) 

As in Sect. 2 it is convenient to use the Feigin-Fuchs free-field realization for the original 
CFT. The space of states (|2.1|) can be realized as 



ncFT = ®p Hp; np = Vp® V-p , (3.11) 

where Vp are the "irreducible Fock spaces" defined in Sect. 2 and V-p are similar spaces 
generated by the "right-moving" chiral Bose field 

m = Q-^p-^J2^ ■ (3-12) 



Each of the spaces Hp forms an irreducible representation of Vir © Vir with the highest 
weight (A(p), A(]3)) and the highest weight vector 



A{p)) =\p)^\ -p) , (3.13) 
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where | p) and | —p) are the vacuum vectors in jFp and J^-p, respectively. The sum in 
is taken over certain set of admitted values of the parameter p which is specific for 
each CFT and determines its Virasoro representation content. In different theories this set 
can be discrete or continuous. 

One remark is in order here. The above perturbed QFT (|3.4|) can be constructed only 
if the original CFT meets certain criteria. For instance, as in CFT the space of local fields 
is isomorphic to the space of states ( |3.11| ), the sum in ( |3.11|) must contain the term with 
2p = 3/3^ — 1 corresponding to the field $1,3 (and its conformal descendants). Moreover, 
according to the fusion rules which follow from the null- vector equations ( p.2|) the operator 
$1^3 acts as 

$1,3 : Hp ^ Hp-p2 (BHp® np+^2 , (3.14) 



and therefore to admit the action of this operator the space (|3.11| ) must contain the whole 
sum 

Ti-p = ®nT^p+nf3'2 (3.15) 

together with each Tip included therein. For particular values of p (corresponding to so 
called "completely degenerate" representations of Vir) this sum truncates and contains in 
fact only finitely many terms (see [10[)- In ciny case we will assume that the space ( |3.11| ) 
has the structure 

nCFT^®anp^. (3.16) 

Obvious example of CFT with this structure is "minimal CFT" Mr/r' with = ^ [|10[ in 
which case p in (|3.11| ) run over finitely many values 2pi^k =1 — k 0^ (/ = 1, 2, r — 1; k = 
1,2, ...,r' — 1). Another example is c = 1 CFT of uncompactified free Bose field where p 
in ( p.ll| ) takes continuous values. 

The perturbation term in (|3.4|) is relevant (i.e. Ai^3 < 1) for Q < (3^ < 1 and so the 
perturbation theory in /x'^ is superrenormalizable (moreover, for < 1/2 this perturbation 
theory does not have any short-distance divergences at all). Therefore in the finite-size 
system with R < 00 where possible infrared divergences are also well contained there are all 
reasons to assume that the space of states TLpcFT of the perturbed theory ( |3.4|) coincides 
with the space of states ( pTTI ) of the original CFT 



TLPCFT — 'HcFT 
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(3.17) 



Of course the highest weight states | A{p)) (|3.13| ) in Hcft are no longer stationary states 
of the perturbed Hamiltonian ( ^.101 ) ; finding the eigenstates and eigenvalues of H in ( |3.17| ) 
is exactly the problem we want to address here. 

In Sect. 2 we have defined the operator valued matrices Lj(A) for the "left" chiral 
sector of the CFT; these matrices satisfy the Yang-Baxter relation (|2.17| ). Similar way one 
can define the "right" chiral counterparts of these matrices, 



L,(A) 



Pexp A / ( : e-^^'^("") : Ej+ : 



Fj) dz 



-inPHi 



(3.18) 



where (j){z) is the right chiral Bose field (|3.12| ), Ej, Fj and Hj denote the same {2j + 1) x 
(2j + 1) matrices as in (|2.13|) , and like in ( |2.13|) V indicates the path ordering of both 
operators and matrices along the integration path. The elements of the matrices ( |3.18| ) 
act in the space JFp defined similarly to ( |2.15|) in terms of the "right" Fock spaces JFp. It 
is possible to check that the operator matrices (|3.18|) also satisfy the Yang-Baxter relation 
in the form 

■/ 

.7 .7 



E E [^^(Wa [L,'(A')]i [R,r{X/X' 

C=-j C'=-j' 



B B' 
CC 



(3.19) 



y: e [^..'(a/v)]^;; [l,,(a')]2; [l,(a)], 

C=-j C'=-j' 

where Rjj'{X) are exactly the same /^-matrices as in (|2.17|) . Note that the i?-matrices 
enter (|3.19| ) different way as compared to ( |2.17| ), the difference being traced down to the 
difference in the commutation relations 



{x),(j){x')] = -[(j){x),(j){x')] = ef^ 



2i '\ R 

where e{x) is the quasi-periodic continuation of the sign function, 

e(x) = 2n+ 1, n<x<n-|-l:nGZ. 



(3.20) 



(3.21) 



Now we can introduce the following "full transfer-matrices" 



T,(^|A)=tr^^.[L,(A)L,(^/A)J = ^ [L,(A)]^ [L,(^/A)] 

A,B=-j 



(3.22) 



where ^ is a parameter. It is easy to check that (|3.22|) defines the operator which acts in 
the space (|3.15|) , i.e. 



T,(^|A) : 



Hp Hp . 



(3.23) 
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It follows from (|2.17| ), (|3.19|) and the fact that the i?-matrices satisfy the "unitarity con- 
ditions" 

E E [R^rW]TA' [RjA^-')\cc' = 5% , (3.24) 

C=-j C'=-j' 

that these operators form the commuting family, i.e. 

[T,(^|A),T,.(//|A')] = (3.25) 



for all values of j,]'. The following simple properties of the operators (|3.22|) are immedi- 
ately established 

1. For |U = the operator reduces to the operator Tj defined in Sect. 2, 

Tj(0|A) = Tj(A) ®I , (3.26) 

i.e. it acts as Tj (A) in the "left" spaces Vp and as the identity operator I in the "right" 
spaces Vp in (|3.11j )i. Similarly, 

T,(/i|/i/A)|^^Q = I®T,(A) , (3.27) 

where I is the identity in the "left" spaces Vp and = tr^. [Lj{X) e'^""^^^ ] act in the 

spaces Vp as the "right" counterparts of the operators (|2.18| ) . 

2. Like Tj(A) the operators Tj(/x|A) are the single-valued functions of A^. 



3. As follows directly from the definition (|3.22|) the operators Tj{^\X) inherit their 



large A asymptotics from those of the operators Tj(A), 

logTj(;u|A) ~ Kj R/{2tx) A^+^ a ^ cx) , (3.28) 

where 

sinfyrvf) 

K = A-^!^ K (3.29) 

^ sin(7r^/2) ^ ' 

and K is given by ( p.23|) . Similarly, the asymptotics of these operators at A ^ are 



controlled by the second factor in (|3.22|) and so 



logTj(;u|A) ~ Kj R/{2n) {fi/X)^+^, A ^ . (3.30) 



^ Obviously, the zero mode operator P, when acting on the space ( p.ll| ), can be identified 
with -P. 
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In particular, it follows from ( p.30| ) that Tj{fj,\X) has an essential singularity at = 0. 

It is clear that for generic value of the parameter /i the operators Tj(/x|A) do not 
commute with the CFT local IM ( |2.4[ ). Instead, it is natural to expect that the operators 
Tj(/i|A) play the same role in perturbed theory ( |3^ ) as Tj(A) did in the CFT, namely, the 
asymptotic expansions of Tj{fx\X) near the essential singularities at A^ = and A^ = oo 
generate the local IM (|3.8| ) of the perturbed theory (|3.4|) . Note that the spaces (|3.15P of 



the perturbed theory admit the action of the perturbed local IM (|3.8| ). The parameter /x 
in ( ^.221 ) is expected to be related to the coupling parameter // in (|3.4|). 
Concerning the operators Tj we can formulate now our basic 

Conjecture 2 

a) The operators Tj{p\X) are single-valued functions of \^ , regular everywhere in the 
complex \^ plane except at = and A^ = oo where they have essential singularities. 

b) The asymptotic expansions of the operators Tj{n\X) near A^ = oo and A^ = are 
controlled by the perturbed local IM ( \3.^ ); in particular, for T{fx\X) = Ti/2(/u|A) one has 



logT(^IA) = KR/{27r) X'+^ - ^ X^'-^-^^'+^^ hn-i , 

(3.31) 

logT(^|A) = «:i?/(27r) (^/X)'+^ -J^ Cn {li/xf hn-i 



where k and Cn are exactly the same constants as in ^2.23( ) and the parameter n is related 
to jl in p.^ j as 

r4(i-/32) / r(3/3^)r(/32) y ^ 

^ " 7r(l-2/32)(3/32-l) Vr(l-3/32)r(l-/32)J ^^"^^^ 



In particular, the "transfer-matrices" Tj commute with all local IM ^3.8^ ) 



[T,(;u|A),l2n-i] = [T,(;u|A),l2n_i] = . (3.33) 



The relation (|3.32| ) appears very natural in view of nearly obvious relation of the 
operators T(|u|A) to the sine-Gordon model (which we nonetheless discuss briefly in Sect. 6 
below); we have arrived at (|3.32| ) by combining ( |2.23| ), ( |3.28|) - (|3.30| ) with known relations 
(6.3), (6.4) and (6.6). 
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It is possible to show that the operators ( p.22|) satisfy exactly the same "fusion rela- 
tions" as the operators Tj, namely 

TJ{^i\q-^X)TJ{^i\q--^X)=I + T^+,/2{^i\X)T^_,/2i^i\X) . (3.34) 

Also, if is a root of unity the fusion relations ( |3.34|) truncate exactly the same way 
as it happens in the unperturbed theory [|l| (see Sect. 4 below). This leads to closed 
functional equations for the eigenvalues of the operators Tj{fx\X) in the space HpcFT- 
Solving these equations one can, in principle, determine those eigenvalues and whence 
find the eigenvalues of all the local IM ( p.8|) in the perturbed theory (|3.4| ), including the 
eigenvalues of the Hamiltonian ( ^.lUI ). In Sect. 5 we will show how this approach works in 



the SLYM (O). 



4. Fussion Truncation and Kliimper-Pearce Equations in CFT 

The properties of the operators Tj described Sect. 2. hold for arbitrary positive (3'^. 
As was discussed in |jl|] considerable simplifications occur if is a rational number (in 
this case q defined in (|2.14| ) is a root of unity and associated value of c corresponds to a 



"minimal CFT" [0). In these cases the "fusion relation" ( p..2| ) truncates and reduces to a 
functional relation involving finitely many operators Tj. This "truncated" fusion relation 
leads then to the closed system of functional equations for the eigenvalues Tj{X) of these 
operators. 

Let g be a complex number of the form 



= ±1 , and q^ ^ ±1 for any integer < n < , (4.1) 



where N > 2 some integer. For such values of q the A^ + 1 dimensional representation 
of Ug{sl{2)) with the spin j = N/2 becomes reducible, while all the finite dimensional 
representation with j < N/2 remain irreducible. More precisely, this A^ + 1 dimensional 
representation reduces to a direct sum of the A^ — 1 dimensional representation with j = 
N/2 — 1 and two one-dimensional representations i. These representations bring separate 
contributions to the trace ( p.l8|) and therefore 

T^^/2(A) = 2cos(27rArp) + T^^/2-l(A) , (4.2) 



^ The latter are special one-dimensional representations p± of Uq(^sl{2)^ with q^^ = 1 for 
which p±{E) = p±{F) = and p±{H) = ±N. 
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where we used the notation P for P to emphasis its operator nature. Due to ( [4 .21 ) the 
fusion relations (|1.2| ) close within operators Tj(A) with j = 0, 1/2, 1, N /2 — 1/2. 

This truncated fusion relation leads to closed functional equations for the eigenvalues 
of the operators Tj(A); j = 0, 1/2, 1, N/2 — 1/2. Let Tj{X) be the eigenvalues of these 
operators associated with some common eigenstate in Vp. The functional equations take 
most convenient form if one makes use of the rapidity variable 

e = log{X'+^R/i2n)) , (4.3) 

and introduces the functions M 



YAO) = 1 (A) T^.+ 1 (A), j = i, 1, f , f - 1 ; 



(4.4) 



Then the functional equations can be written as 



{o + '-f) Y, [e-'-f) = {i + Y^_.{e)) (1 + Y^^.{9)) , , = 1, 1, . . . , f - 1, 



YN_Ae + — ]YN_Ae — ^ 



2 / -~'\ 2 

= {l + YN_s{e)) {l + e'^'P^Y{e)) (l + e-2'^^P^F(6')) 



(4.5) 



F(« + !|£)F(«-!^i) = (i + r„_.(«)) 



where as before p denotes the eigenvalue of the operator P. Note that the system 
coincides with the functional form of the TBA equations of I? at type PB| . 

Further simplifications occur when p takes special values. Suppose that the operators 
Tj act on the space Vp with 

£ + 1 

p = —A- , (4.6) 
where £ > is an integer such that 2p ^ nfP + m for any integers n and m. Then one can 



show ||21|| that in addition to (|4.2| ) the operators Tj satisfy 



T^_,_i(A) = (-l)^T,(A), for = 0, i, 1, f - 1; 

(4.7) 

Tw_i(A) = . 

2 2 ^ ' 

This is exactly the case for the "minimal CFT" M.2/2n+z with n = 1, 2, where 

27ri 

q = e2r.+3 
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and the space of states has the form 



2/2ra + 3 



(4. 



with 



Pk 



2/c + l 



/c = 0, 



2(2?i + 3)' ' 2(2n + 3) 

In this case the extra relations (|4.7|) aUow one to bring ( [4 .51) to yet simpler form 

Y,{e + Z7r/(2n + 1)) ¥,{9 - Z7r/(2n + 1)) = (l + Y^_Ae)) (l + F.+ i (^)), 



(4.9) 



1 



(4.10) 



Again, (|4.10|) coincides with the functional form of the TBA equations for the perturbed 
"minimal CFT" M 

2/2n+3 |]23- The above truncation was discussed in 
In this paper we concentrate attention on the simplest of the above "minimal CFT" , 
the model M.2/5- This CFT describes the criticality associated with the Lee- Yang edge 
singularity and therefore it is often referred to as the Lee- Yang CFT. The central 
charge is 

22 

c(-M2/5) = -y (4.11) 
and the space ( [4.8| ) contains only two components 



where po = 1/10 and pi = 3/10, the associated conformal dimensions are 

A{po) = -1/5, A(pi) = . 



For this model the functional equation (|4.10| ) takes the form 



Y{e + z7r/3) Y{6 - z7r/3) = 1 + Y{e) 



(4.12) 



(4.13) 



(4.14) 



where according to ( [1.4| ) the function Y{6) = Yi{6) = Yi{6) simply coincides with corre- 
sponding eigenvalue Ti/2{X) of the operator T(A) = Ti/2(-^) in the space (|4.12| ), 



F(^)=Ti/2(A) , e' = XsR/{2n). 



(4.15) 



The analytic properties of the eigenvalues Ti/2{X) was discussed in Sect. 2. They are 
entire functions of the variable 



t = XU — 

2tv 



R \ 5 
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with the asymptotic behavior 



logTi/2(A) ~ K te, as t^oo, (4.16) 

where 

"' '^'/T (r(3A))' = 28.29877111.... 



as it foUows from ( |2.23D with /3 = 2/5. Correspondingly, we are interested here in the 



solutions Y{9) of ( [4.14| ) which are regular everywhere in the complex ^-plane, satisfy the 



periodicity condition B 

Y{9+^i7:)=Y{e) , (4.17) 

and have the asymptotic 

y(^) ~exp(Ke^), as ^ ^ +oo . (4.18) 

Note that since the eigenvalues Ti/2{X) are regular at t = the associated functions Y{9) 
are bounded at ^ ^ — oo; in fact in this limit they approach the constant values 

Y(n\ _^ f 2cos(27rpo) = (1 + "\/5)/2 for the eigenvalues in Hp^ ; (A ^Q\ 
\ 2cos(27rpi) = (1 — y/5)/2 for the eigenvalues in Hp^ . 



as it follows from ( P.26q| ) 



The functional equation identical to ( 4.14| ) arises in the solvable lattice hard hexagon 



model [H. It was studied by Klumper and Pearce [^, who developed the method of 
transforming the functional equations to nonlinear integral equations, the knowledge about 
analytic properties of the solutions being key ingredient in this approach. Although in the 
field theory the analytic properties of the functions Y{9) are somewhat different from those 
considered in 0, the method of Kliimper and Pearce can be easily adopted to our case. 
The derivation in is based on the following 

Lemma. // a function f{9) is regular and bounded in the strip '^m9 E (— 7r/3, 7r/3) 
and satisfies the relation 

f{9 + m/3) + f{9 - ztt/S) - f{9) = g{9) (4.20) 



^ In fact, this periodicity condition can be derived from ( [4.14| ), see [22|. 
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with some g{6) then 

/oo 
de' ^{e-e')g{e') , (4.21) 
-00 

where 

^ ^ « . (4.22, 

TT sinh(36') 



Note that 



^[e) = ^de\ogS{e) , (4.23) 



where S{6) is the S-matrix ( p..4| ). If l^(^) was free of zeroes in the strip 53m 6 G (— 7r/3, 7r/3) 
then the function f{6) = logy(^) — k exp(^) would satisfy the conditions of this lemma 
with g{9) = log (l + Y~^{9)), so the integral equation would follow immediately. However 
a typical eigenvalue T{t,p) does have zeroes in the wedge — 27r/5 < argt < 27r/5 (which 
corresponds to the above strip) and so we need to take care of them. 

Let us assume that Y{9) has + 2M zeroes in the strip Qm9 G (— 7r/3, 7r/3), real 
zeroes a = 1, 2, and M complex-conjugated pairs (/3b+z76, (3b—i'yb), b = 1,2,..., M. 
We introduce the following functions 

ao(^)=tanh(^^) , 
cosh(6') — cos(?7) cosh(6') — sin(7r/6 — ?]) (4.24) 
cosh(6') + cos(?7) cosh(6') + sin(7r/6 — ?]) ' 

which satisfy the equations 

ao{9 + in/3) ao{9 - Z7r/3) = 1; ai{9 + Z7r/3, ry) ai{9 - z7r/3, 77) = ai{9, ry) , (4.25) 
and write Y{9) as 

N 

Y{9) = \[ao{9-aa) e<'^ 

a=l 

N M 

= exp(/^ e^) \[a^{9 - aa)X[cT^{9 - ^b.lb) X{9) . 

a=l 6=1 

In view of ( [4.25| ) the function X{9) satisfies the equation 



(4.26) 



N 



X{9 + iTv/3) X-\9) X{9 - iTv/3) = JJ ao{9 - a^) (l + Y-\9)) . (4.27) 



a=l 
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Note that because the cr-factors in ( 4.26| ) absorb all zeroes of Y{9) in the strip 9 G 



(— 7r/3, 7r/3) and the exponential prefactor saturates the asymptotic (|4.18|) the function 



logX{9) is regular and bounded in this strip. Applying the above lemma we obtain the 
integral equation 

M .no N 

49') 



d9' ^9-9') log ( l[ao{9' -aa) + 

0=1 a=l 

(4.28) 



where e{9) is defined in (|4.26|) 



Given a set of real parameters aaif^b^lb the equation ( ^.28| ) uniquely determines the 
function e(^) and hence Y[9). Of course these parameters can not be taken arbitrary. 
With generic values of these parameters the function Y{9) defined through the solution 
to ( |4.28D will have the poles in the domain $5m^ G (tt/S, 27r) mod 57r/3. Indeed, if, for 
instance, Y{[3b + i'^b) = then according to ( [4.14| ) the function Y{9) must exhibit a pole at 
9 = (3b + — 27rz/3, unless, of course, Y{f3b + 27^ — z7r/3) = —1. Therefore the condition 
that Y{9) is entire function leads to the equations 

Y{aa + iTT/3) = Y{aa - itt/S) = -1, a = 1, N ; 

(4.29) 

Y{(3b + ilh - ^7r/3) = Y{^b - + ^7r/3) = -1, b= 1, M . 

which determine the parameters aa, (3b and 7^ in (|4.26|) . The integral equation ( [4.28| ) 
together with the transcendental equations ( [4.29|) are QFT versions of the equations ob- 
tained in the lattice theory in and therefore we will refer to them as Kliimper-Pearce 
equations. 

The Kliimper-Pearce equations ( [4.28| ), ( ^.29| ) enables one to determine (in principle) 
any eigenvalue Y{9) of the operator T(A) in the space ( [4.12| ). In this paper we consider only 
two simplest eigenvalues - the vacuum eigenvalues T(t,po) and T{t,pi) which correspond 
to the vacuum vectors \po > and \pi > in the spaces Vpp and Vp^ in (|4.12|) , respectively. We 
will use the notations Yo{9) = T{t,po) and ^1(6') = T(t,pi) (which are not to be confused 
with the used in (0) - (gT^)) with t = exp{69/5). 

I. Yo{9). This eigenvalue corresponds to the ground state in the space ( |4.12|) . As 
discussed in Sect. 2 all zeroes of T(t, 1/10) lie on the negative part of the real t-axis and 
so all zeroes of Yo{9) are located along the lines Qm9 = 5n/6 (mod 5n/3). In particular, 
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there are no zeroes in the strip QmO G (— 7r/3, 7r/3). Therefore in this case the equation 
( [4.28|) takes the form 

/oo 
de' $(6' - e') log(l + e-^°(^')) , (4.30) 
-oo 

where eo{9) = logYo{9), and the equations ( |4.29D do not appear. The Eq.( ^.3U| ) is exactly 
the massless TBA equation associated with the Lee- Yang CFT The solution of this 
equation was referred to as the "kink solution" in where it was studied numerically; 
the solution is plotted in Fig.l. 

8 r 

6 - 

4 - 
2 - 





^ ^ ^ ^ ^ ^ ^ ^ ^ ^ : 

-4 -J -2 -1 Q 1 2 

norm 

Fig.l. The functions Yo{9) and yi(^) found from numerical solutions of the integral 
equations ( [4.30| ) and (4.37) plotted versus a "normalized" rapidity variable ^norm = 
6 + log K. 




Note that for 6 — oo the function 10(6*) approaches the value 2cos(7rpo) = (v^ + 
l)/2 = 1.6180339. . . in agreement with ( |4.19| ). We have used the numerical solution to 
the equation ([4.30| ) to evaluate few first coefficients in the expansion (|2.26a| ) ; the result is 
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Exact values 


Numerical values 


G^«c(-L)= 1.6180339887... 
^^"^(-L) = 7.01811993... 10^ 
0^''%^) = 1.361347 ... 10^ 
G'^«c(^) = 1.6320... 10^ 
GT^(l^) =(unknown) 


1.618033989. . . 
7.0181198... 10^ 
1.36135... 103 
1.6317... 10^ 
1.41... 10^ 



Table 1. The exact eigenvalues of a few first NIM for the A = — i vaccum state 
in the Lee- Yang CFT (as given by (A. 2) and (A. 9) with a = 2/5 and p = 1/10) and 
numerical values of the same eigenvalues obtained with a polynomial fit of Yo{6) deter- 
mined from the numerical solution of the integral equation ( [4.30 ). The quantity Gq'^'^{p) 
denotes the constant term in the series expansion (|2.26aD 



given in Table 1 where the numerical results are also compared with the exact coefficients 
calculated as the integrals ( p.27|) (see Appendix A). 

As is explained in Sect. 2, the function €o{9) —k exp(6') can be expanded into asymptotic 
series in powers of exp(— ^). The easiest way to obtain this expansion is to note that the 
kernel $ in (|4.30|) expands for 3?e ^ > as 



m 



oo 
n=l 



sm 



27r(n+ 1) 



,(l-2n)6 



(4.31) 



Substituting this into ( [4.30|) one finds 

log yo(^) - ^ - ^ f; sin xn e^^"^")^ , 

n=l 

where the coefficients Xn are calculated as 

Xn= I de e(2--i)^ log(l + e-^o(^)^ 



(4.32) 



(4.33) 



According to ( p.266|) the coefficients in ( [4. 321) are related to the vacuum eigenvalues of the 

local IM l2n-l, 

(4.34) 



ttC. 



2 27r(n + l) 
— sm Xn , 
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where 



^""3 n! UJ r(2i^) ^^^^/^^ 



(4.35) 



The first of the integrals ( |4.33| ), can be evaluated exactly using the well known "dilogarithm 
trick" (see Appendix C) 

. (4.36) 



TT 



15 ^ 



This agrees with the value /J'"'^(— 1/5) = —1/60. We used the numerical solution to 
( [1.30|) to evaluate few further coefficients ( [4.33| ); these results are compared with the exact 
eigenvalues /2^1]^(— 1/5) in Table 2. 



Exact values 



Numerical values 



j_ 

'60 
89 



756000 

211 _ 
8100000 

_ 2160997 
464373000000 

6283403 



1924560000000 



— 1 . 6666666666666 
= 1.17724867724 
= -2.6049382716 
4.6535802 
-3.264851 



..10-2 
.. lO-'^ 
. . 10"^ 
. . 10-^ 
. . 10"*^ 



'I7l 



(unknown) 



-1.6666666666665 . 
1.1772486773, 
-2.604938279 . 
4.653584. 
-3.26488. 
3.5501. 



10-2 
10-4 
10-^ 
10-^ 
10-*^ 
10-^ 



Table 2. The exact vacuum eigenvalues of a few first (non- vanishing) LIM given 
by (B.1)-(B.8) with c = —22/5 and A = —1/5 and numerical values for the same LIM 
obtained from ( [4.34 ) with the numerical solution of the integral equation ( [4.30| ). 



II. Yi{9). According to our discussion in Sect. 2 all zeroes of the function T(t, 3/10) are 
real, just one of them being positive. Correspondingly, the function Yi{9) = T(t, 3/10) has 
one real zero which we denote a and infinitely many zeroes located along the lines Q^m 6 = 
57r/6 (mod 57r/3). Therefore in this case the function ei{9) = logYi(6') — logo"o(6' — a) 
solves the equation 

e^{9) = Ke^ + f d9' $(6* - 9') log (tanh(3(6'' - a)/4) + e-^^^^'^) , (4.37) 
and the parameter a satisfies 



I e 



El (a+i7r/3) 



-I e 



ei (a — i7r/3) 



1 . 



(4.38) 
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The last relation can be written as 



r- 3 cosli(2(^-a)) / 

v^Ke" + --/- d9 ^ log ( tanh (3(^ - a)/4) + e" 

TT sinli(3(6' — ajj V 



J-oo sinli(3(6' — a 
-7r(l + 4iV) , 



(4.39) 



where N is some integer and -;f_^ denotes the principal value of the singular integral. We 
have used ( |4.37| ) and explicit form of the kernel $(6') (\i.22\i to obtain (|4.39|) . We studied 
the equations ( [1.37| ) and ([4 .391 ) numerically. Our analysis suggests that these equations 
admit real solution only if in ( ^.39| ) is zero, in which case 



a + log K = 0.49577315... . 



The corresponding function 



Yi{e) =ao{0-a)e 



(4.40) 



(4.41) 



is plotted in Fig.l. In this case the asymptotic value of ^1(6*) at 6* —00 is 2 cos(27rpi) ~ 
(1 - v^)/2 = -0.61803398. . .. Few first coefficients of the expansion ( p.26a|) calculated 
from this numerical solution together with their exact values (when available) are presented 
in Table 3. 



Exact values 



r^vac( 3 \ 


-0.6180339887... 


-0.618033989. 




r~<vac( 3 \ 







3.6.. 


10"^ 


r~<vac( 3 \ 


2.94659 . . . 


102 


2.9467. 


..102 


r^vacf 3 N 
<^3 \W - 


5.80714... 


10^ 


5.80705 . 


. . 10^ 


r^vacf 3 \ 


6.2826. . . 


10^ 


6.281. 


. . 10^ 


r~fvacl 3 \ 

^5 - 


4.729. . . 


10^ 


4.75. 


. . 10^ 


r^vacl 3 N 

^6 lioJ - 


2.72. . . 


10*^ 


2.6. 


. . 10^ 


r^vacl 3 \ 

^1 - 


(unknown) 


1.7. 


..10^ 



Numerical values 



Table 3. The exact eigenvalues of a few first NIM for the A = vacuum state 
in the Lee- Yang CFT (as given by (A. 2), (A. 7), (A.8), (A.IO) with a = 2/5 and p = 
3/10) and numerical values of the same NIM obtained with a polynomial fit of Y\{Q^ 
determined from the numerical solution of the integral equation ( 4.37| ). The quantity 



Go"'^(p) denotes the constant term in power series expansion (2.26a) 
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To obtain the asymptotic expansion for log Yi(6') similar to ( [4 .321 ) one could try again 
to substitute the expansion ( ^.31] ) into ( [4.37| ). However in this case the log factor in the 
integrand in ( [4.37| ) does not decay sufficiently fast as 6*' ^ oo and so the integrals appearing 
after this substitution typically diverge. To handle this problem it is convenient to rewrite 
the equation (|4.37| ) as 



log \Yi{9)\ = K log {S{e -a + in/3) S{9-a + 2in/3)) + 

/oo 
de' ^{e-e') iog( | i + Yf\e') \ ) . 
-oo 

Now using (f4.31D we obtain the asymptotic series expansion 



logyi(6') = K e > sm 

TV ^ — ^ 3 



Xn e 



(l-2n)0 



where 

Xn 



2tv 27r(n+l) 
sin 



271-1 



/oo 
d^e(2-i)^log(|l + yfi(^)|) 
-oo 



(4.42) 



(4.43) 



(4.44) 



Numerical values 



Exact values 



1.83332. .. 10-1 
4.5106. . . 10-4 
-5.56789... 10-5 
7.2631... 10-6 
-4.717. . . 10-6 
4.67. . . 10-6 



/i(0) = 
^5(0) = 
/7(0) = 

^11 (0) 
/i3(0) 
-^17(0) = (unknown) 



11 

60 

341 
756000 

451 
8100000 

^ 3373117 
464373000000 

^ 825517 
174960000000 



1.833333, 
4.51058, 
-5.567901 , 
= 7.26380, 
= -4.7183, 



. 10-1 
. 10-4 
. 10-5 
. 10-6 
. 10-6 



Table 4. The exact vacuum eigenvalues of a few first (non- vanishing) LIM given by 
(B.1)-(B.8) with c = —22/5 and A = and numerical values for the same LIM obtained 



from (4.46) with the numerical solution of the integral equation (4.37) 



Like in the previous example, for n = 1 the integral in ( [4.44|) can be calculated 
analytically using appropriate modification of the "dilogarithm trick" (see Appendix C), 

XI = -^^-'. (4.45) 
15 
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The vacuum eigenvalues /2ri-i(0) °f local IM are expressed through Xn as 



2 27r(n + l) 
— sm Xn , 



(4.46) 



where Cn are given by ( [4 .351 ), so that (|4.45[ ) agrees exactly with /J'"^(0) = 11/60. The 
values of few further local IM obtained from the above numerical solution are compared 
with the exact values in Table 4. 



5. Excited States in Scaling Lee- Yang Model 

In this section we apply the approach outlined in Sect. 4 to the finite size energy 
spectrum of the SLYM (|1.3| ) . The QFT ( |1.3| ) is particular case of ( |3.4| ) and so it possesses 
infinitely many local IM ( |3.8| ). In fact, due to the strong degeneracy of the representations 
entering ( [1.12|) in this case some of these IM vanish. Only the IM l2n-i with 2n — 1 7^ 
(mod 3) are nonzero. As is mentioned in the Introduction the finite-size ground state 
energy Eo{R) of this model was studied by TBA in Q]. In this section we show that the 
approach based on the operators ( p.22|) and the fusion relations ( p.34|) makes it possible 



to study the excited state energies in the finite size system as well. 

As was explained above, the perturbed theory inherits the space of states from the 
CFT M2/5, that is it has the form 

where po = 1/10 and pi = 3/10. Consider the operator T{fx\X) = Ti/2(a*I-^): as defined in 
( p.22|) . This operator invariantly act in the space ( |5.1|) and when specialized to this space 
it satisfies the truncated fusion relation 

T(^|QU)T(^|Q-iA) =I + T(^|A) , (5.2) 



and therefore the eigenvalues T(//|A) of T(//|A) in the space ( |5.1| ) satisfy exactly the same 
functional equation 

T{^i\qh)T{^\q-^X) = l+T{fi\X) (5.3) 

as the eigenvalues of the operator T in the space ( [1.12| ). However, in the perturbed case 
(/Lt 7^ 0) we are interested in the solutions to ( p.3| ) which have significantly different analytic 
properties as the functions of A^. As was explained in Sect. 3, the eigenvalues T(/x|A) have 
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essential singularities both at = oo and = 0, and they exhibit the asymptotic 
behavior dictated by ( p.28|) and (|3.30|) near these singularities. 



Like in the Sect. 4 it is convenient here to introduce the rapidity 

e = ^ log(AV^) . (5.4) 

Note that this definition differs from ( [4.15|) by a constant shift log(/ut R/2'k). Also, we 
will use the dimensionless variable 



r = m 



R= {KiJ^/n)R (5.5) 



instead of R. Here m is the mass related to the coupling parameter /} (|1.3|) as in ( |1.5|) . 
The eigenvalues T{n\X) can be considered as the functions of 9 and r. Like in ( [4.15| ) we 
denote 

Y{r\e) =T{fi\X); r = KfiiR/7r, = . (5.6) 



In this notations the functional equation ( |5.3| ) takes the form ( [4.14| ), i.e. 



Y{r\e + z7r/3) Y{r\e - iiv/S) = 1 + Y{r\e) . (5.7) 

We are interested now in the solutions which are analytic everywhere in the complex 9- 
plane, satisfy the periodicity condition ( ^.17| ) and have the asymptotic 

logy(rl^) ~ r , 9 ±oo . (5.8) 



Like in the conformal case discussed in Sect. 3 the solution to ( ^.7|) with the asymptotic 
conditions ( p.8| ) is completely determined by the pattern of zeroes of Y{r\9) in the strip 
— 7r/3 < Qm9 < n/3. The same transformations as in Sect. 3 lead to the integral equation 



M 



e{r\9) =r cosh ^ + ^ log ai [9 - Pb, 7b) + 

(5-9) 

d9' ^9 - 9') log ( n MO' - tta) + e-^("l'') 

a=l 



where 

n 

e{r\9) = logy(rl^) - logao(^ - a„) (5.10) 



a=l 
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and we have assumed that there are + 2M zeroes in the physical strip, real zeroes 
aa{r) and M complex-conjugated pairs (3b{r) ±z76(r); here we have shown explicitly that 
the positions of these zeroes depend on the parameter r H. And again, Uke in Sect.4 these 
positions are determined by the equations 

y(r|a„ + z7r/3) = y(r|a„ -z7r/3) = -1, a = l,...,iV; (5.11) 

Y{r\(3i, + ilb-i^/^) = Y{r\(3t-nb + iTr/3) = -l, 6=1,...,M. (5.12) 

Before studying particular examples of the eigenvalues Y{r\9) let us show how the mass 
spectrum and factorizable S-matrix of this QFT can be recovered from these equations in 
the limit r oo. Simple estimate shows that for r >> 1 the first term in the r.h.s. of 
Eq.( ^.9|) dominates everywhere in the strip Qm9 G {—n/3 + e , 7r/3 — e) except in the small 
disks of the size e around the zeroes and Pb^i'Jb, where e is some number exponentially 
small in r. In this domain the first term in the r.h.s of (|5.7D can be neglected and hence 
the function Y(r\9) can be written as 



Y{r\e) ^YlS{e - Pa - in /2)exp{r cosh 9) , 

a 

r » 1; \9 - f3a±in/6\> e; Qm9 e {-n/3 + e,n/3 - e) 



(5.13) 



where >S'(^) is the S-matrix defined by (|1.4[ ) and /3a are some real numbers. In writing 
( p.l3| ) we have taken into account the asymptotic conditions ( [5.8| ) and the fact that Y{r\9) 
must be real at real 9 i. The form ( p. 13 ) shows in particular that for r sufficiently large 



there are no real zeroes in (|5.9| ) while the complex zeroes have the form 

(3b{r)±i{n/6-Ub{r)) , (5.14) 

where Pa{f) are real and z/a(r) are exponentially small in r. The equation ( |5.12| ) then takes 
the form 

exp(ir smhPb) II S{Pb- Pa) = I, 6=1, 2,..., M, (5.15) 



We will see later that when r change some complex-conjugate zeroes can colhde and then 
become a pair of the real zeroes, or vice versa, so that in different domains of r the numbers A'^ 
and M may vary, while the total number of zeroes, of course, remains the same. 



The factors of the form S{6 + m/6) S{9 + 5i7r/6) also satisfy (5.7) with the first term in the 
r.h.s. in this equation omitted, and they agree with the real analyticity condition. More detailed 
analysis shows however that these factors would be incompatible with ( ^.12| ). 
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where we have taken into account that 5'(0) = — 1. This equations shows that for r >> 1 
the eigenstate associated with the eigenvalue ( ^.131) is interpreted as M-particle state, the 



parameters (3b = (3b{r) being the rapidities of the particles, so that the equations ( |5.15| ) 
have the meaning of the Bethe Ansatz equations for M scalar particles with the S-matrix 
( p. .41) . We should stress here that in the field theory the Bethe Ansatz equations of 
the form ( |5.15| ) make sense only in the limit r —>■ oo. For finite r these equations acquire 
corrections due to the vacuum polarization. 

Here we consider in details only two simplest eigenvalues of the operator T of the 
Lee- Yang model. 

I. The simplest eigenvalue is of course the ground state eigenvalue, i.e. the eigenvalue 
associated with the lowest eigenstate of the Hamiltonian ( p.lO| ) in the space (|5.1| ). We 
denote here this state as | \E'o) and use the notation Yo{r\9) for the corresponding eigenvalue 
of T. For the ground state eigenvalue Yo(^l^) there are no zeroes in the strip 9 G 
(— 7r/3, 7r/3) at all and the equation ( |5.9| ) takes the form 



/oo 
de' ^{6 - e') log(l + e-^o('^l^')) , (5.16) 
-oo 

where eo(r|^) = loglo(^l^)- This is exactly the TBA equation for the Lee- Yang model 
obtained in 0]; here we have arrived at (|5.16| ) by entirely different route, completely 



bypassing any reference to the on-shell solution of the model. The equation ( 5.16 ) was 
studied in very details in [|[ and we can add very little to that analysis. Let us just briefly 
summarize the results of [^. For large r the flrst term in the r.h.s. of ( |5.16D dominates 



and so eo(r|^) ~ rcosh^; r >> 1. In the opposite limit r ^ the function yo(^l^) 
becomes nearly constant equal to in the wide range of 6*, logr < 9 < — logr, while 

for 6* > — logr this function coincides with the "kink solution" Yq(6 + \og (^)) considered 
in Sect. 4; similarly, for 6 < logr, YQ{r\6) Yq{^ — 6 — log (^))- In view of ( |3.26| ) and 
( p.27|) this is just the manifestation of the fact that for r — > the ground state | \E'o) 
approaches the conformal vacuum state | A(]3o)) ( |3.13| ) in (|5.1|) . According to (|3.31| ) the 



ground state eigenvalues of the local IM 

l2n-l |^o)=41l(i?) I^'O) 

i2n-i |vE'o)=/S-i(^) 1*0) 
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(5.17) 



can be obtained from eo{r\9) as the coefficients of the asymptotic expansion 

oo 

eo{r\e) ^ reV2 - ^^'^'-'^y^ e^'-'^^' I^tiiR) , ^ - +00 , 

n=l 
00 

earn ^ re-V2 -f^C^ ^5(i-2n)/6 ^(2n-i). ^(o)_^(^) ^ ^ ^ -00 . 



(5.18) 



It foUows from ( |5.16| ) that 



/r)(i?) = -||m^i?-£ |^d^e^log(l + e--H^)) , 

= ^ f— sin^!^i!?±ll r e(--)^ log(l + e--('-|^); 



(5.19) 



where n = 2,3, ... . The similar relations hold for the "right" local IM Then, 
for the ground-state energy 

H I ^'o) = Eq{R) I ^'o) , 

we have 

Eo{R) = ^(i?) + i[''\R) = -— m^R - — / dO cosh6' log (l + g-^o^'^l^)) . (5.20) 

12 27r J— 00 

Obviously it is the first term in ( |5.20| ) that dominates at large R and therefore it gives the 
bulk vacuum energy of SLYM while the second term represents the finite-size corrections. 
Finally, the leading R ^ asymptotic 

Eo(R) ~ (5.21) 

can be deduced from ( |5.20| ) in full agreement with expected form Eo{R) = — (7r/6i?)(c — 
24A(po)). 

II. Let us consider now the first excited state of the SLYM. We will denote this state 
and associated eigenvalue of T as | and Yi{r\9), respectively. At ^ 00 this state 
is interpreted as one-particle state with the particle momentum (and rapidity) equal zero. 
Therefore, according to our discussion above, for large r the eigenvalue Yi{r\9) has two 
zeroes in the physical strip located at 



9 = ±i 7(r) = ±i (tt/G — z^(r)) , 
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(5.22) 



where z^(r) — as r — oo with exponential accuracy. Correspondingly, in this case the 



integral equation (|5.9| ) for the function ei{r\9) = logyi(r|^) takes the form 

/oo 
d6' $(6* - 6') log(l + e-^i(''l^')) , (5.23) 
-CXD 

where ai{9,r]) is given in ( [4.24| ) and the function z/(r) which describes the deviation of 
zeroes from their limiting positions at ±in/6 is determined from the equation (|5.12| ) 

gei(r| ±i(7i-/6+i/)) ^ _-|^ _ (5.24) 

In the limit r >> 1 the term with the integral in ( |5.23| ) is exponentially small and with 
this accuracy one has 

cosh y + cos(g- — z^J coshc/ + smz/ 
where u = u{r) is derived from ( [5.24| ) 

z/(r) = V3 e" " + O (re"^") . (5.26) 

For r ~ 1 the equations ( [5.23| ), (|5.24| ) can be solved numerically. The resulting function 
7(r) is given in Table 5 for 2.6 < r < 8.0 and the typical form of the function Yi{r\9) in 
this range of r is shown in Fig. 2. 




r 


7(r) 


AE{r)/m 


(TCSA) 


2.600000 


0.133055 


0.238050 


0.23805 


2.800000 


0.256131 


0.197403 


0.19741 


3.000000 


0.322401 


0.164321 


0.16432 


3.200000 


0.367132 


0.137217 


0.13722 


3.400000 


0.399638 


0.114891 


0.11489 


3.600000 


0.424188 


0.096412 


0.096411 


3.800000 


0.443176 


0.081057 


0.081055 


4.000000 


0.458106 


0.068256 


0.068251 


4.200000 


0.469989 


0.057552 


0.057552 


4.400000 


0.479533 


0.048579 


0.048583 


4.600000 


0.487254 


0.041042 


0.041040 


4.800000 


0.493538 


0.034700 


0.034699 


5.000000 


0.498676 


0.029356 


0.029352 


5.200000 


0.502895 


0.024846 


0.024847 


5.400000 


0.506370 


0.021037 


0.021034 


5.600000 


0.509241 


0.017816 


0.017812 


5.800000 


0.511619 


0.015092 


0.015091 


6.000000 


0.513592 


0.012785 


0.012782 


6.200000 


0.515233 


0.010832 


0.010837 


6.400000 


0.516598 


0.009178 


0.009177 


6.600000 


0.517737 


0.007775 


0.0077767 


6.800000 


0.518688 


0.006587 


0.0065876 


7.000000 


0.519482 


0.005581 


0.0055816 


7.200000 


0.520147 


0.004727 


0.0047264 


7.400000 


0.520703 


0.004003 


0.0040045 


7.600000 


0.521169 


0.003389 


0.0033970 


7.800000 


0.521559 


0.002870 


0.0028753 


8.000000 


0.521886 


0.002429 


0.0024364 



Table 5. The data obtained from the numerical solution of ( ^.23| ) and (5.24). The 
third column contains the excess energy AE\(r) defined by (5.36) (given in units of 
the mass m) while the forth column contains the TCSA data for the same quantity 
computed in JH 



As r decreases the zeroes (|5.22D get closer together and at 



ro ~ 2.53576 



(5.27) 
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they collide on the real ^-axis. For r < vq the function Yi{r\9) has two real zeroes 

6 = ± a{r) , (5.28) 

and one can represent it in the form: 



Fi {r\e)= ao {e - a{r)) ao {9 + a{r)) e 



r < To 



' tanh ( 



3a(r) 



2 4 
instead of and (^M )- 



r 


a{r) 


AE{r)/m 


(TCSA) 


(CPT) 


1.010"^ 


10.39926 


23037.48 




23037.348 


1.010-3 


8.096681 


2302.837 




2302.834 


1.010-2 


5.794092 


229.3856 




229.3849 


0.100 


3.49150 


22.0527 




22.0527 


0.500 


1.88203 


3.67974 




3.67977 


1.000 


1.18800 


1.44738 




1.44736 


2.000 


0.464431 


0.430871 


0.43088 


0.430197 


2.400 


0.206153 


0.288398 


0.28839 


0.286762 


2.500 


0.103278 


0.261863 


*0.26189 


0.259921 


2.510 


0.0875353 


0.259338 


*0.25939 


0.257393 


2.520 


0.0684264 


0.256826 


*0.25691 


0.254891 


2.525 


0.0565558 


0.255570 


*0.25569 


0.253650 


2.530 


0.0414288 


0.254311 


*0.25447 


0.252416 


2.535 


0.0150627 


0.253036 


*0.25325 


0.251188 



(5.29) 



where the function ao{9) is given by ( [4.24]) . Now we should solve the equation 

/oo 
de' ^e-6')log{ao{e'~a{r))ao{e' + a{r))+e-^'^^'^\^'^) , (5.30) 
-oo 

together with 



Table 6. The data obtained from the numerical solution of (|5.30| ), ( |5.31| ). The 
third column contains the excess energy A£^i(r) defined by (5.36) (given in units of the 
mass m) while the forth column contains available TCSA data for the same quantity 
computed in |f3[|(the numbers marked with an asterisk obtained by the qubic spline 
interpolation of the resutls of |p^). The fifth column contains CPT results for the same 
quantity obtained from (5.37). 
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0.06 



0.0-1 - 



o.o: 



0.00 



-0.02 - 



-0.04 I I 

2.30 2.40 2.50 2.60 2.70 

Fig. 3. The values of a'^(r) and — 7'^(r) for r ~ ro and the hnear fit ( |5.32D (dashed 

hne). 

The function a{r) obtained by numerical integration of these equations is given in 
Table 6 for 10~^ < r < vq. Note, that the location of zeroes of Yi{r\9) for r ~ tq both for 
r > ro and r < ro is described by the same law (see Fig. 3) 

^ ^ ^ A(ro - r), r ~ ro, (5.32) 



with the numerical value of constant A = 0.297893 
For r ^ the zeroes ( p.28|) depart to ±oo as 



a(r)~-log(;^)+a, r^O, (5.33) 



where the constant a is the same as in ( [4.40p . In this limit the function Yi{r\9) approaches 



the constant 2cos(27rpi) = (1 — V5)/2 in the wide range of 6, logr < 6 < — logr, while 
yi(r|^) ~ yi(^ + log (^)) for ^> - logr and Yi{r\e) ~ Fi ( - - log (^)) for 0< logr. 
Here ^1(6*) is the "kink solution" discussed in Sect. 4. A typical form of the function Yi{r\9) 
for r < ro is shown in Fig. 2. 
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J 



- ^ 




(a) (b) 

Fig. 4. Diagramms representing two leading contributions to the large r asymptotics 
of A£'i(r). Diagramms (a) and (b) coreresponds to the third and fourth terms in (5.35), 
respectively. 

The eigenvalues of the local IM \in-\ associated with the excited state | can be 
calculated from the asymptotic expansion of ei(r|6') for r > ro and ei(r|6') for r < ro 
similar to ( p.l8|) . In particular, for the energy Ei(R), 



this gives 



Ei{R) ^ mr + 2m sm( - + u{r)) - — / d0 coshO log(l + Yr^{r\0)), r > tq , 

12 '6 ' 2tv 

Ei{R) — mr + -\/3m cosh a (r) / d6 cosh. 6 \og\l -\-Y^^{r\6)\, r < tq , 

12 27r 

(5.34) 

Substituting ( |5.25D and ( [5.26| ) into the first of the above formulae one obtains a few first 
terms of the large r = mR expansion. 



Ei(R) = mr + m + 3me~~^- 

^ ' 12 



/OO 
de cosh 6* S{e + m/2) e-'-'^^^he ^ oire-^"^) 
-OO 



(5.35) 



where S{9) is given by (|1.4|) . Here the first term is the bulk vacuum energy, the second 
term makes it manifest that for R ^ oo the state | describes one particle with zero 
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-6 
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0.1 
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Fig. 5. The dashed Hne represents the CPT correction terms (the second and third 
terms in (5.37) with negated signs given in units of the mass m) and the points represent 



the values for ll7r/(15r) — Ei{R)/m obtained from numerical integration in (|5.3 



momentum, while the third and the fourth terms are the corrections coming from the 
diagrams in Fig. 4. These terms in ( |5.35| ) are in fuU agreement with the results of Q , 
and [|7|. 

In the opposite limit R ^ Q the state | approaches the CFT highest weight state 
I A(pi)) (|3.13|) in (|5.1|). The later has the dimensions (A, A) = (0,0) and corresponds 
to the identity operator in this CFT. Indeed the calculation in Appendix C shows that 
Ei{R) ~ ll7r/(15i?), as — s> 0. In Tables 5 and 6 the function Ei{R) obtained by 
numerical integration in ( |5.34| ) is compared with the first excited state energy obtained in 
by Truncated Conformal Space Method. To facilitate this comparison we present the 
data for an "excess energy" 



\/3 

AEi(r) = Ei{R) + ^mr 



m, 



(5.36) 



where the bulk and mass terms are subtracted. Also, in Table 6 we compare our numerical 
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results with the first terms of the short-distance expansion 



Ei{R) = m(^ + e2r^+e3r^+---); (5.37) 
lor 



e.--(l + V5)*|i^™(r(2/3)r(5/6))-; 



TTir i-(l/5) 

29 ^JL ^3, 
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obtained by Conformal Perturbation Theory from ( |1.3| ) (see Appendix D). In Fig. 5 we plot 
numerical results for Ei{R)/m — ll7r/(15r) against the correction terms in 



6. Discussion 



In Sect. 3 we have constructed the commuting family of T-operators for perturbed 
CFT (|3^). As is well known the IQFT can be obtained by so called "quantum group 
restriction" of the sine-Gordon model . The sine-Gordon model is described by 

the action 



AsG = 



1 



167r 



{daifY - 2^* cos{(3ip) 



d X 



(6.1) 



where {p{z,z) is a scalar field and (3 and fi are parameters, the latter one carrying the 
dimension 



H ~ [length'\ 



2/3 



(6.2) 



This QFT is integrable, in particular it possesses infinite set of local IM of the form (3.8) 
where now T2k+2, 02fc and T2k+2, 02fc are certain local fields of the sine-Gordon model 
(and fi is replaced by fl). The on-shell solution of ( |0|) contains two "topologically charged" 
particles — soliton and antisoliton — and a number of the neutral particles with the masses 



2Msin(7rjO, j 



2'"^''" ^ 2^ 



(6.3) 



^ Note, that for r < 0.1 the correction terms in ( 5.37 ) are quite small (less than 10^^ of the 
magnitude of the leading 1/r term) and are beyond the precision of our numerical calculation. 

In fact the value and even the dimension of the parameter jl depends on the precise way 
one defines the composite field cos(/3(/?) in ( |6.lD . Here we assume that this field is canonocally 
normalized with respect to its short-distance asymptotic, i.e. cos{Pip){z,z) cos{f3(p){z' , z') ^ 



1 \z 

2 1^ 



[length] 



{z,z) {z',z'). In this normalization the field cos{f3ip) carries the dimension 



and hence the /i in ( |6.lD has the dimension (6.2). 
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where ^ is defined in ( p.24|) and M is the sohton mass which in turn is related to the 
parameter Ji in (|6.1|)as [|lj 



M 



2 r(|) 



i±i 

2 



(6.4) 



The factorizable S-matrix of the sine-Gordon model can be found in IRC 



In the case of infinite space the model ( |6.1[ ) exhibits the symmetry with respect to the 
quantum group Uq{SL{2)) 



29|] where 



(6.5) 



Namely, the soliton and the antisoliton transform as the two-dimensional representation of 
this quantum group and the local IM ( p.8| ) and the S-matrix commute with the generators 
E,F and H of the associated quantum algebra Uq{sl{2)). Naturally, the Hilbert space 
Hoo of the infinite-space sine-Gordon model containes the subspace H^^^^^ of the states 



anihilated by the above generators. The remarkable fact discovered in [28] , |2!]] is that 



this subspace H^^^^^ can be interpreted as complete space of states of certain local QFT. 
It turnes out that this "restricted sine-Gordon model" coinsides with the perturbed CFT 
(pl^), the parameter fi in (|3.4| ) being related to the constant ft in ( |6.1|) as 



(l-2/?2)(3/32_l) 



TT 



r3(/j2)r(i-3/3^ 
r3(i-/32)r(3/32 



(6.6) 



Although the action of the generators E, F, H can be defined also in appropriate sectors of 
the finite-size sine-Gordon model with the spatial coordinate compactified on a circle, the 
quantum group symmetry does not survive this compactificationS. Nontheless the notion 
of singlet states still exists, and the subspace '}{^^^3iet ^ ^^-^ g^^j^j^g interpretation 



as the space of states of the perturbed CFT 



/jisinglet 



in finite size geometry 



'HpCFT 



(6.7) 



The generators E,F,H can be defined only in particular "good" sectors of the space of 
states Ti.R of the finite-size system, with appropriately chosen "twists" (i.e. the eigenvalues of the 
operators dxf and dyip. This is similar to the supersymmetric QFT where in the finite-size 
case the SUSY generators exist only in the Ramond sector. However unlike the SUSY theories in 
general case of q the generators of Uq{sl{2)) change the "twists" and do not invariantly act in the 
"good" sectors. 
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It is clear from the structure of the operators Tj{^\X) defined in ( |2.13| ), (|3.18| ) and ( |3.22| ) 
that their action extends to the fuU sine-Gordon space Tin where they also satisfy the 
commutativity conditions (|3.25|) and all further relations in Sect. 3 including ( p.31|) and 
( p.34|) hold in the full space. In other words the operators Tj{^\X) can be considered as 
the QFT "transfer-matrices" of the full unrestricted sine-Gordon theory. From this point 
of view the fields (p{z) and (p{z) admit very simple interpretation in terms of the sine- 
Gordon field, namely cp and (p are just the values of the field <f>{z, z) on the "left" and 
"right" components of the light cone, respectively, i.e. 



-|(^(0,i? 



(6.8) 



here we find it more convenient to think in terms of the Minkowski space-time, where 



X — t, 



(6.9) 



are the real light-cone coordinates. Then according to ( |3.22| ) the operator Tj{fx\X) is 
obtained by integrating the Lax fiat connection associated with sine-Gordon model (see 
e.g. [0 ) around the space-time cilinder, the integration being done in two stages — first 
one integrates along the "left" light cone from z = to z = R with z kept equal to 0, and 
then the integration along the "right" light cone from z = to z = R is performed,^! as is 
shown in Fig. 6. 




Fig. 6. The integration path corresponding to the form ( ^.221 ) of the sine-Gordon 
"transfer matrix". 



This choice of the integration path makes calculations simpler. As the connection is flat one 
could have chosen any other closed contour which wraps around the cihnder. 
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It is possible to check that if takes special values such that (^4.1|) is satisfied the 
truncation relation ( [4 .21) holds in the full sine-Gordon space of states (the relation ( [4 .71 ) 
is specific to the restricted theory) and so an appropriate version of the method outlined 
in the Sectins 4 and 5 above can be applied to the sine-Gordon model at these values 
of However for generic values of the approach based on the QFT version of the 
Baxter's Q-operator is expected to be more powerfuU. We have shown in how to 
construct the Q operator in CFT, as the trace similar to ( |2.18| ) but this time taken not 
over a finite-dimensional reppresentation of Uq{sl{2)) but over certain infinite-dimensional 
representation of so-called (/-oscillator algebra. Clearly, the operator Q(|u|A) for the sine- 
Gordon model can also be defined by the formula similar to (|3.22| ). As we observed in 
||20[] the asymptotic expansion of logQ(A) in CFT contains the contributions of the "dual 
nonlocal IM" in addition to the local IM. It is easy to see that the "dual nonlocal IM" are 
closely related to the nonlocal IM discovered in sine-Gordon theory in . Therefore the 
operators Q(/u|A) and associated Baxter's equations are expected to be a powerfull tools 
in calculationg the spectra of both local and nonlocal IM in sine-Gordon theory lil. Our 
progress in this direction is under way p^ . 



Our construction for the T-operators can be extended to the cases when higher rank 
quantum algebras replace Uq{sl{2)). In the case of CFT the operators T introduced in 



was recently generalized to the case of g-deformed twisted Kac-Moody algebra A2 in ||33|| . 



The "massive" versions of these operators similar to (|3.22|) would apply to another class 
of IQFT — the CFT perturbed with the operator $1^2 — which includes in particular the 
T — Tc Ising QFT with nonzero magnetic field |TI . 



In fact, the relation ( |3.32| ) in our Conjecture 2 in Sect. 3 is obtained by combining (2.23), 
dU) and (U). 
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Appendix A. 

In this Appendix we use the notation 



a = . (A.l) 



While the simplest of the integrals (p.27| ) G\°''^{p) is easily evaluated 



^^"^(p) = / dx dyi 2 sin { ) j 2 cos [2p {x - y - tt)) 

47r2r(l-2a) 



47r sin (7r(a - 2p)) [ {1 - t)-^" t"'^^-^ dt 
Jo 



r(l -a-2p)r(l -a + 2p) ' 

(A.2) 

no explicit expression for these integrals with /c > 1 is known for generic p and a. The 
integral G2"''^{p) can be simplified as 

Grip) = (2.r cos(2.p) ^'(^ - r(« - 2P) r(a + 2p) 



r(l -a-2p)r(l -a + 2p) 

o -p2/i ^ / sin7r(a- 2p) r{2p) sin7r(a + 2p) r{-2p) 

— 27iri 1 — a ; — -, -r(a,p) + - ^ — ; -r(a,—p) 

^ ' ^ {a-2p)r{l-2a-2p) ^ {a + 2p) T{1 - 2a + 2p) ^' ^' 

(A.3) 



where the function F{a,p) is given by the series 



F(a,p) = r(l +20) ± „'t:!CrL /"(«.P) . 



with 



n! r(3a — 2p + n) 



(x — p 

/o(a,p) = l, /i(a,p) = 8a— — , (A.5) 

1 + a — 2p 
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and the rest of the coefficients /n(fl,p) with n > 1 determined recurrently through the 
relation 

An{a,p) fn+i{a,p) - Bn{a,p) fn{a,p) + Cn{a,p) fn-i{a,p) = , (A.6) 

where 

An{a,p) = {l + n + a-2p){l + n- 4p); 

Br,{a, p) = 2n^ + I2n^ {a - p) + 2n (8p^ + 12a^ - 4ap - 7a) + 8a (4p^ + a - iap - p); 

Cn{a,p) = 71 (n + 4a — 1) (n + 3a — 2p — 1) (n + 4a — 4p — 1) . 

This expression is very efficient for numerical evaluation of G2"''^{p) with high precision. 
Another convenient expression for this integral which also generalizes to higher is 
given in 

In particular case 2p = 1 — a the integrals G2"''^{p) and G'^°''^{p) can be evaluated 
analytically 

rvac^-^^ 2vr^ r3(l-a)r(3-4a) 
2^2^ l-2a r2(2-2a)r(2-3a) ' ^ ' 

G^ac(l^^ ^ r3(l-a)r2(l-2a)r(f-3a) 
3^2^ l-2a r(2-4a)r(2-3a)r3(f -a) ' ^ 



Finally, note that when g is a root of unity the functional equations (^) or ( ^4.1U| ) 



impose certain algebraic relations between the nonlocal IM, which allow one to express 
some higher nonlocal IM through the lower ones. In particular, using the functional ( [4.15| ) 
one can show that 

5-^5, 



(A.9) 



for the p = 1/10, A = — 1/5 vaccuum state in the Lee- Yang CFT and 



'-^4 lloJ - I 10. 



Grm = arm (A.io) 



5 - 

To" 

nvaci 3 \ _ ^ "t" ( r^vac( 3 ^ ( nvac/ 3 \ 

^6 KW - -^Q \W) ^ — iToJ 

for the p = 3/10, A = state. 
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Appendix B. 

Here we present eigenvalues of the local IM l2n-i = I7 2, 8) with R — 27t on the 
highest Virasoro vector with the central charge c and conformal dimension A 



ir(A) = A - ^ , (B.l) 

,~c(A) = A._(£±^A+£(5£±M, (B.2) 

^ ^ ^ 12 2880 ' ^ ^ 

/^-^(A) = A3 - i£±il + + + A - + , (B.3) 

^ ^ ^ 8 576 290304 ' ^ ^ 

jvac(^s ^ ^4 _ (c + 6) ^3 (15c^ + 194c + 568) _ (c+2) (c + 10) (3c + 28) ^ 
^ ^ ^ 6 1440 10368 

c (3c + 46) (25c2 + 426c + 1400) 
^ 24883200 ' ' 

^ ^5 _ M£±8) ^4 ^ (c + 8)(5c + 46) 

^ ^ ^ 24 288 

(35c3 + 990c2 + 9048c + 23488) ^2 
48384 

(c + 2) (175c3 + 7134c2 + 96168c + 392000) ^ 
11612160 

_ c (5c + 22) (11c + 232) {7c^ + 27Ac + 1960) 
3065610240 

rvac(.. ^ . 6 _ £+10 ^5 , 15c^ + 322c + 1808 ^4 
^ ' 4 576 

_ 105c3 + 3700c2 + 44612c + 165984 

72576 " 

525c4 + 27908c3 + 548508c2 + 4248784c + 10147200 , o 

H 

11612160 

(c + 2) (315c4 + 24604c3 + 676548c2 + 7298480c + 25872000) 



+ 



(B.5) 



,3 



418037760 

c (13c + 350) (11025c^ + 1160780c3 + 25741404c2 + 198779728c + 470870400) 

27389834035200 ' 

(B.6) 
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r-r ^ A 7 _ 7(c+12) 6 7(15c^ + 386c + 2680) , 
^3 ^ ^ 24 2880 

(105c3 + 4430c2 + 66264c + 319552) ^4 
41472 

(525c4 + 33364c3 + 8181720^ + 8304848c + 27924096) ^3 
^ 4976640 

(3465c5 + 331192c4 + 120126680^ + 1909592960^ + 1291772608c + 2856307200) ^2 

1313832960 

+ (c + 2) (I21275c^ + 18838640c^ + 9229340360^ + 188032083520^ + 166162628800c 
+ 517957440000) A/3310859059200 

- c (5c + 164) (3465c^ + 1026934c^ + 390094760^ + 5680476560^ + 3512182240c 

+ 7399392000) /79460617420800 , (B.7) 

r-f Al = A« - ^^+^^^ A^ + 7(c^ + 30c + 248) , 

(7c3 + 344c2 + 6140c + 36304) ^5 
1728 

(35c4 + 2588c3 + 760200^ + 952528c + 4102528) ^4 
165888 

(231c5 + 25654c4 + 11224760^ + 22259304c2 + 195500000c + 597940480) ^3 

32845824 

+ (21021c^ + 3835642c^ + 236079916c^ + 64389519280^ + 842012499200^ 
+ 500631533696O+ 1026155648000) AV143470559232 

- (0 + 2) (3003c^ + 11378140^ + 890281480^* + 28522628560^ + 441003664640^ 

+ 324626771840O + 897792896000) A/1721646710784 
(5o + 22) (7o + 68) (17o + 658) 
^ 28097274319994880 ^ 

(4290^^ + 4957240^ + 200213880^ + 285527760o + 1457456000) . (B.8) 



+ 



For = —2 the eigenvalues of the aU local IM is known [20|: 



iZ-i(^) = B2n{2p + 1/2) , with A = 2p2 _ 1/8 , (B.9) 



where Bk{x) is the Bernoulli polynomials . 
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One can also find the eigenvalues /2^^i(A) for c = 1/2 and A = 0, 1/2, 1/16: 

^2:-l(0)L=l =-Cn (1-21-2-) 



)l-2n\ -^2n 



4n ' 



/2n-i(l/2)L^. =c. (2^-2- -(1-2^-2-) ^) , (B.IO) 



4\"-i(l/16)L_i=c. 
Here is the Bernoulli numbers and 



(6?i - 3)!! n\ 2"-i 
~ (2n- 1) (4n-3)! 3" ' 

Appendix C. 

In this Appendix we present the calculation of first integrals in ([4 .331 ) and ( |4.44| ) 
(corresponding to the ground and the first excited state in the Lee- Yang CFT) using the 
so-called "dilogarithm trick" . Although for the ground state case such calculations are well 
known in the literature on TEA (see e.g. for the most transparent presentation) we 

include them here for a completeness. 

Introducing the notation 

Lo(^) = log(l + e-^°(^)) , (C.l) 
we can write the first integral in p.33| ) as 

/oo 
de e^Loie) . (C.2) 
-oo 



"OO 

Xi 



Differentiating the integral equation ( [1.30|) with respect to 6, expressing therefrom the 
term and substituting the result into ( |C.2| ) one obtains 

/oo 7*00 7*00 

de deeoie) Lo{e) - de de' ^{e - e') Lo{e) de'Lo{e') . (c.3) 
-oo J — oo J — oo 

Using now the integral equation ( [4.30| ) to simplify the second term in the last formula one 
can rewrite it as 

/oo /"OO 
de {deeo{e)Lo{e)-eo{e)deLo{e))-K / dee'Lo{e). (C.4) 
-oo J — oo 
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Thus, 



de {deeo{e)Lo{e)-eo{e)deLo{e)) 



A change to a new integration variable 



x{9) = (l + e^°(^)) ' 



(C.5) 



(C.6) 



brings this this integral to the form 



Xi K = —Li (x(oo)) + Li oo)) 



(C.7) 



where Li (x) is the Rogers dilogarithm function 



and 



x(— oo) = ( — ^ , x(oo) = 



as it follows from ( ^.ISp and ( [4.19| ). Here and below we will need the following special 
values0 of Li (x) 



Li(l) = — , Li, 
^ ' 6 ' V 2 



— , Li 
10' 



~ 15 



(C.9) 



Using this in ( |C.7| ) one arrives to the result ( [4.36|) given in the main text. 

Let us now turn to the calculation of the first integral in ( [4.44| ) . It can be conveniently 
rewritten as 

dOe^Liie), (C.IO) 



Xi 



where 



Li(^) =log (ao(^-a) + e-^i(^)) 



(C.ll) 



with ao{0) defined in ( [4.24| ). Proceeding as above we can bring ( |C.1C1| ) to the form similar 
to dCj) 

1 

(C.12) 



^ I -oo 



14 



See, e.g. [37| for a comprehensive review on the dilogarithm function. 
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(C.13) 



(C.14) 



(C.15) 



Defining new functions 

m = ei(^) + log \ao{e -a)\= log , 

L{e) = Liie) - log |ao(^ - a)| = log |1 + Y-\e)\ , 
one can further rewrite (|C.12|) as 

X^K = I+^J de (dem L{e) - m deL{e)) , 

where / is determined as the principal value of the singular integral 

de [Li{e)+ei{e))de\og\ao{e-a)\ . 

-oo 

The integral in ( |C.14| ) is calculated similarly to that in (|C.5| ). Splitting the interval of 
integration into two parts, (—00, a) and (a, 00) and changing to a new integration variable 

e^^^\ for 9 < a 

one can bring ( |C.14|) to the form 

Xi K = / + Li (1) + Li ( x{-oo) ) , (C.17) 

where Li (x) is the dilogarithm function ( |C.8|) and 

V5 - 1 



x{9) 



x{—oo) 



(C.18) 



as it follows from ( [4.19| ). The remaining integral (|C.15| ) is calculated by using the integral 
equation ( |4.37| ) and the condition ( f4.39| ) determining the root a. In fact, using ( [4.37| ) first 
write ( |C.15|) in the form 

/oo POD 
l dOde' Li{e) {^e- 6') + 5{e-e')) de' log \ao{e' -a)\ , (C.19) 
-00 ^ — 00 

where we have used the elementary integral 



d6 e^deloglaoie - a) 



(C.20) 



The integral over 9' in ( |C.19| ) is also elementary and can be easily calculated. After that 
the r.h.s. of (|C.19|) becomes (up to a simple factor) identical to the l.h.s. of ( [4 .391) 

y_oo smh3(6'-a) 
CoUecting together (|C.17| ), ( |C.21|) , (|C]|) one obtains 



(C.21) 



/22 \ TT^ 

( — +24A^) — 
V 5 / 6 



XiK = -[ — +-z^i\ ) — , (C.22) 

which for = precisely gives (|4.45| ). The above calculations can easily be generalized 
for the solutions of ( |4.28| ) and ( f4.29| ) with arbitrary number of (real and complex) zeroes. 
Similar calculations for in lattice models can be found in HI. 
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Appendix D. 

One obtains Conformal Perturbation Theory (CPT) by expanding in the interaction 
parameter ji^ in ( p.4|) . The CPT for the ground-state energy Eq{R) of ( |1.3| ) was discussed 
in [Q. It can be easily modified for the excited-state energies. The excited -state energy 
Ei{R) in Sect. 5 expands as 

^ 27r \2fc(A-i)+i 

X 



^ ^ 6R ^ k\ \ R J 

k=2 



fD.ll 



where (i^z = dx dy for z = x+iy, c = —22/5 is the central charge of CFT A = — 1/5 

is the conformal dimension of the field $ in ( |1.3|) and (. . ■)'cft (connected) 
correlation functions of the this CFT 111. The first terms in ( p.l|) contain 

mZi,ZiMz2,Z2)rCFT = \Z1-Z2\-'^ , (D.2) 
($(^l,^l)$(^2,^2)$(^3,^3))cFT = C$$<I. \zi - Zsl"'^ 1^2 - ^231"'^ 1^3 " ^1 1" • (D.3) 

Here ^ 

c$$$ = iT^ r3(i/5)r(2/5) 

is the triple-^ OPE structure constant []T^. The most important difference in ( p.l| ) as 
compared with the similar formula for Eo{R) in is that the integrals in ( p.l| ) diverge 
as {zi,Zi) — > 0, oo and therefore the expression (|D.1|) can not be taken literally. This 



divergence is easily traced down to the fact that | is not a ground state. As is well 
known in ordinary quantum mechanical perturbation theory, in doing calculations for the 
excited states in each order one has to take care of the orthogonality between the state 
under consideration and all the lower-energy states. This condition brings into ( p.l| ) 
additional terms which exactly provide the subtractions of the above divergences. So 



in (p.l|) appropriate subtractions of these divergences are implied. Equivalent way to 



define the integrals in ( p.l| ) (which is similar to "dimensional regularization" of Feynman 



15 



For generic excited-state energy En{R) the expansion similar to ( p.l ) would contain the 



insertion Ki(0, Q)Vn{oo, oo) in the correlation functions, where Vn{z, z) is the CFT fields associated 
with the CFT hmit of the excited state | under consideration. For the state | '^\)r=q the 
associated field is the identity operator and therefore this insertion is ignored in ( p.l]) . 
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diagrams) is to evaluate the integrals in ( p.l| ) for A < 1/2 where the integrals converge 
and then analytically continue to A = —1/5, see e.g. [^. With this prescription the 
integrals in the first two terms of the series in (|D.1|) can be evaluated explicitly, 



r^(A)r(i-2A) 
' r2(i-A)r(2A) 



(D.4) 



d zi d Z2 \zi 



2A-2 



-2A 



TT 



, ^ll "|1 

r^(f)r(i-^) 



Z2\ 



-2A 



Zl - Z2\ 



-2A 



(D.5) 



4 r3fi 



3A^ 



2 / V 2 



and one arrives at ( |5.37| ). 
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